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Motivated by ongoing measurements at JILA, we calculate the recoil-free spectra of dipolar in- 
teracting fermions, for example ultracold heteronuclear molecules, in a one-dimensional lattice of 
two-dimensional pancakes, spectroscopically probing transitions between different internal (e.g., ro- 
tational) states. We additionally incorporate p-wave interactions and losses, which are important 
for reactive molecules such as KRb. Moreover, we consider other sources of spectral broadening: 
interaction-induced quasiparticle lifetimes and the different polarizabilities of the different rotational 
states used for the spectroscopy. Although our main focus is molecules, some of the calculations are 
also useful for optical lattice atomic clocks. For example, understanding the p-wave shifts between 
identical fermions and small dipolar interactions coming from the excited clock state are necessary 
to reach future precision goals. Finally, we consider the spectra in a deep 3D lattice and show how 
they give a great deal of information about static correlation functions, including all the moments 
of the density correlations between nearby sites. The range of correlations measurable depends on 
spectroscopic resolution and the dipole moment. 
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I. INTRODUCTION 

Recoil-free spectroscopy is a powerful, pervasive probe 
of ultracold atomic systems that includes radio-frequency 
and microwave spectroscopy, as well as Doppler-free op- 
tical or Raman spectroscopy. Interatomic interaction ef- 
fects on the spectra give information about the phase dia- 
gram and behavior of many-body systems. For example, 
these effects on recoil-free spectra were used to first de- 
tect Bose-Einstein condensation in dilute spin-polarized 
atomic hydrogen pQ, locate and probe the Mott insula- 
tor/superfluid quantum phase transition of bosons in op- 
tical lattices [2H5J , and study Cooper pair binding [7J [5] , 
polaron quasiparticle residue [5] , and pseudogap behavior 
of ultracold fermions across the BEC/BCS crossover [TDT - 
I12j . In precision measurements, e.g. atomic clocks, un- 
derstanding the effects of interactions is crucial because 
they can limit [13] or enhance [T3J the achieved accuracy. 

Here we calculate the recoil-free spectra of dipolar 
molecules in a one-dimensional lattice of two-dimensional 
"pancakes" and in a three-dimensional lattice. Such sys- 
tems have been realized at JILA [L5HI9] . As we demon- 
strate, recoil-free spectroscopy gives information com- 
plementary to density profile or time-of-flight measure- 
ments. Additional probes are especially important for 
molecules, where direct imaging is extremely difficult due 
to the absence of cycling transitions [20] . To be concrete, 
we consider driving the transition between two rotational 
states of the molecule, but the same calculations apply to 
transitions between vibrational, hyperfine, and electronic 
states when the transitions are measured with techniques 
imparting negligible recoil momentum. 
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There is great excitement about ultracold molecules: 
they offer additional internal degrees of freedom (rota- 
tional, vibrational) compared to atoms, and the long- 
range character of their interactions can lead to interest- 
ing many-body physics |21H24j . Examples are quantum 
liquid crystals [2"5rl2"7] . Wigner crystals [3SJ[2D], exotic su- 
pcrfluids [30 32j , supersolids [55J , topological phases [5il - 
157] . and the physics resulting from quantum phase tran- 
sitions among this plethora of phases [351 EH] ■ We expect 
recoil-free spectroscopy to play a major role in observing 
and characterizing these systems. 

The outline of our paper is as follows. Sec. |III A gives 
general expressions, valid even for strongly correlated 
systems, relating the average spectral shifts to simple 
static correlation functions. Sec. IIII Bl evaluates these ex- 
pressions explicitly for weakly interacting homogeneous 
systems. Sec. |III C| applies these expressions to obtain 
the trapped system's spectra, while Sec. D| shows ex- 
ample results. These calculations are directly relevant 
for the ongoing experiments at JILA [TSlfTD] and should 
be useful in a number of related ongoing efforts to make 
ultracold molecules, e.g. [3D] . 

We quantitatively calculate the spectral shifts and, in 
Sec. |IV B| lineshapes due to the dipolar interactions, p- 
wave interactions - including p-wave losses that are im- 
portant for reactive molecules such as KRb - and the dif- 
fering polarizabilities of the different molecular rotational 
states that may be used in the spectroscopy. We also 
qualitatively consider and bound other effects, such as 
higher lattice band occupation and the validity of treat- 
ing the atoms in the cold collision regime, i.e. treating 
interactions via a pseudopotential. 

To briefly motivate the utility of these spectra, we 
point out that measuring them in the dilute thermal gas 
allows one to ensure that these various non-trivial con- 
tributions — many of which are just beginning to be 



explored in the context of cold atoms — behave as ex- 
pected. Furthermore, deviations from the simple dilute 
thermal gas values allow one to diagnose more interest- 
ing behavior as the temperature is lowered, the simplest 
example probably being the degree of quantum degener- 
acy of the gas. The onset of strongly correlated quantum 
phases will also have signatures. 

As another application of recoil-free spectra to dipolar 
interacting systems, Sec. [V] considers a deep 3D lattice. 
We show that one can measure the joint probability dis- 
tribution P(l,m,7i2, ■ ■ •) °f a particle occupying a site, 
while having n\ molecules in the nearest neighbor sites, 
ri2 in the next-nearest neighbor sites, and so on. This 
includes all moments (n^n^ 2 ■ ■ ■) for arbitrary otj. The 
spatial range of correlations that can be measured is set 
by the spectral resolution and interaction strength, and 
is typically a few sites for the molecules presently consid- 
ered. This provides a powerful way to characterize the 
system, measuring a much more complete set of corre- 
lations and with greater spatial range than other tech- 
niques, such as modulation spectroscopy [4"T] . 

While our quantitative estimates of parameters are for 
dipolar molecules, many of the calculations are also rele- 
vant for Rydberg atoms and alkaline earth atomic clock 
experiments. As one example, Rydberg atoms have long 
range dipolar interactions. As another, p-wave interac- 
tions can be important in optical atomic clocks, as re- 
cently revealed by Ramsey spectroscopy of Yb lattice 
clocks with the atoms trapped in ID arrays of 2D pan- 
cakes or 2D arrays of ID tubes [42] . Moreover, as ex- 
periments in these clocks expand their capabilities to 
approach expected 10 -18 accuracy and mHz frequency 
resolution, small dipolar interactions between the atoms 
become relevant. We expect our calculations may pro- 
vide insight into these experiments. However, we mention 
that while our theory assumes that few atoms are trans- 
ferred to the excited state (i.e., linear response), many 
clock interrogation techniques transfer a substantial frac- 
tion of the atoms [4"31 133] . 



II. THEORETICAL AND EXPERIMENTAL 
BACKGROUND 

We are interested in heteronuclear molecules in a one 
dimensional optical lattice of two dimensional pancakes 
in the presence of an electric field. We assume that the 
energy scale set by this field is large compared to the 
hyperfine splitting; we explain this in more detail be- 
low. Fi g. [TJ illustrates this system and defines notation. 
Section [V|considers the effects of adding a lattice to the 
transverse directions. This system displays a hierarchy of 
energy scales, which dictates the physical behavior. Wc 
begin by describing this hierarchy. 

Figure [2] shows the lowest energy rotational levels of 
molecule in a DC electric field. This structure is well- 
described by the rigid rotor Hamiltonian in an electric 
field, H = BN 2 — d ■ E with N the angular momentum 




FIG. 1: Schematic of the experimental geometry consid- 
ered. Heteronuclear dipolar molecules are trapped in a one- 
dimensional lattice of two-dimensional pancakes, with the lat- 
tice periodic in the z direction. The angle 6 is the angle be- 
tween the inter-molecule separation r — r' and the molecular 
dipole (applied electric field axis), here taken to be along the 
z-axis. 



operator, B the rotational constant, and d the dipole 
moment operator. This single-molecule physics sets the 
largest energy scales in the problem. In the absence of an 
electric field, there is the usual rotational level structure: 
a ground state, a 3-fold degenerate spin-1 first excited 
state manifold, a 5-fold degenerate spin-2 second excited 
state manifold, etc. These are labeled by their total an- 
gular momentum M and z-axis projection m z . In the 
presence of an electric field E the rotational symmetry 
is broken. Here we take E to be along the z direction, 
perpendicular to the pancakes, as is used in the ongo- 
ing KRb experiments to suppress the reactive losses [15] . 
The electric field mixes states with different M but the 
same m z . 

Figure [2] shows the resulting level structure for a non- 
zero value of E, where we label by \M, m z ) the state that 
is adiabatically connected to the E = state with angu- 
lar momentum M and z-projection m z . An electric field 
induces a level splitting v in the previously degenerate 
excited state manifold, and this energy scale is large com- 
pared to the other energy scales that will be of interest. 
As a consequence, when we treat the interactions, lattice, 
and trap, it is quantitatively accurate restrict oneself to 
processes occurring in the resonant subspace where we 
neglect transitions that require energy changes on the 
order of v. 

We now present the Hamiltonian describing the period 
lattice, the harmonic confining potential, and intermolec- 
ular interactions for molecules in an electric field. We 
consider two molecular states, choosing |0, 0) and |1,0) 
to be concrete, but the vast majority of our calculation is 
general and straightforwardly extended to arbitrary in- 
ternal states. We will denote these states as \a) with 
a = 1,2 henceforth. The hyperfine interaction energy 
is the next largest energy scale after the rotational and 
electric field, and can couple the a = 1,2 and other ro- 
tational states, shown in Fig. [2j However, even a small 
electric field - one that generates dipole moments more 
than an order of magnitude smaller than the permanent 



3 



|2, -2) 



|2.-1) 



1 2,0} 



|2,+1) 



|2,+2) 



11,0} 



...Iv 



\M = 0. 



0} 



FIG. 2: Lowest rotational energy levels of a single molecule 
in an electric field. These are described by the rigid rotor 
Hamiltonian in an electric field H = BN 2 — doE with N the 
angular momentum operator, B the rotational constant, and 
do the dipole moment along z. The label \M, m z ) indicates the 
state that is adiabatically connected to the E — state with 
angular momentum M, and angular momentum z-component 
projection m z . However, the E-Re\d mixes states with the 
same m z but with different M, so M is not a good quantum 
numbers. The energy splitting v increases with increasing 
electric field and is a large energy scale even for small E: 
interactions, the trapping potential, and the kinetic energy 
cannot mix in states that are off-resonant by this energy scale. 



dipole moment - suffices to generate level splittings ~ v 
that are larger than the hyperfine coupling. Thus for 
even such small electric fields, if the initial populations 
are confined to these two rotational states, we can ne- 
glect the other states, which are then far off-resonant 
compared to the hyperfine energy. 

The lattice energy scale is the largest after the ro- 
tational, electric field, and hyperfine structure of the 
molecules, all set by single molecule physics. The lat- 
tice separates the energy levels into bands, and because 
the other energy scales - trap, interactions, and tempera- 
ture - are small relative to the band splitting for the deep 
lattices of present interest, we can thus project the Hamil- 
tonian to the lowest lattice band for states a = 1, 2. We 
neglect terms that are exponentially smaller in the lattice 
depth (e.g. tunneling terms). The molecular system is 
described by the Hamiltonian 



H — H + H sw + i/dipolc 

with the non-interacting part 



V 2 

2m 



-Strap (1) 

( 2 ) 



where i/j a (i,p) and ijr a (i,p) are fermionic annihila- 
tion and creation operators acting on atoms in rota- 
tional state a at site i in the lowest band at trans- 
verse two-dimensional position p. They satisfy the 



anti-commutation relationships {ip a (i, p), ipg(J, p')} — 
SijSapS^p — p'). Here m is the molecular mass and fi 
is the chemical potential. 

The projected s-wave interaction Hamiltonian is 

H sw = J2^Y~ J d P^(hP)i>l(hP)^p(hP)^ a (hP) 

(3) 

with g a p — (47ra Q( g jm) J dz \w a (z)\ 2 \wp(z)\ 2 (we set h = 
ks = 1 throughout), where a a p is the s-wave scattering 
length for scattering between molecules in states a and 
/3, and w a (z) is the lowest band Wannier function for 
particles in state a. This depends on a because the lat- 
tice depths felt by each species need not be equal since 
the molecular polarizabilities differ. The interaction gn 
is undefined for identical fermions, which don't scatter in 
the s-wave channel, and so we can set it to zero, but g\2 is 
nonzero. When the Wannier functions for the two states 
are not identical, the three-dimensional p-wave interac- 
tion projected onto the lowest band gives a renormaliza- 
tion of the two-dimensional 1-2 s-wave interaction (not 
included above) in addition to the contribution coming 
from the three-dimensional s-wave scattering, because 
the particles effectively become distinguishable. How- 
ever, the s-wave interaction will turn out to be irrelevant 
for the recoil-free spectra in the linear response regime 
with all molecules initially in the ground internal state, 
as a consequence of the Pauli exclusion principle. Hence 
we can ignore the s-wave contact interaction. 

The dipolar interaction Hamiltonian in 3D is [45, 4G 



V 



(3d) _ 



drdr' 



1 - 3 cos 2 



J"""" 4ire Q \r 
42^iPi(r)<4(r')Pi(r')P2(>") + H - c - 



mm^m m 



,(4) 



where 9 is the angle between the dipole orientation di- 
rection and r — r', d mn = (m|d(°)|n) where d^ is 
the m z — spherical component of the dipole oper- 
ator in the basis aligned with the electric field, and 
ip a (r) and ^(r) are fermionic annihilation and cre- 
ation operators satisfying the anticommutation relations 
{<p a (r), if a(r')} = 5 a p5{r — r'). The reason why only the 
m = spherical component of the dipole operator is re- 
tained is that even a small electric field drives the m = ±1 
levels off-resonance, as shown in Fig. [2} To be sufficiently 
off-resonant compared to the other energy scales - inter- 
actions, lattice, and hyperfine coupling (to allow one to 
neglect the nuclear degrees of freedom) - requires only a 
small field, where the induced dipole moment is still an 
order of magnitude or more smaller than the permanent 
dipole moment. The same reason justifies why only the 
shown matrix elements of the dipole interaction are kept 
- the other matrix elements leave the resonant subspace. 
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Note that the permanent dipole moments of these levels 
are non- vanishing, in contrast to the zero-field rotational 
levels, because the electric field induces a dipole moment 
to the eigenstates. 

The projected dipolar interaction Hamiltonian is 



-^dipole = ^ X! [ d >P d P ,V <x0( i -j>P-p') 

ij,a0 

x ipl(i, p)^(j, ff)il> p (j, f/)rf> a (i, p) 
+ \Y1 j dpdp'V si (i -hp- P') 

ij 

x (4(i, p)4(j, p')^(j, p> 2 (*, P) + Hx.)(5) 

with 

V a p{i- j,p- p') = Idzdz'- 



d aa dpf}(l — 3 cos 2 



47re r 



(p-p'f + iz-zf 



X \w a (z - Zi)\ 2 \wf}(z' - Zj)f 



3/2 



(6) 



dzdz' - 



^12^21(1 — 3 cos 2 



and 

V s {(i-j,p- p') = , - 

47T6Q [( P - p') 2 + (Z - z'Yf' 2 

x w\(z — Zi)w 2 (z — Zj)w* 2 (z! — Zj)w 1 (z' — Zj).{7) 

We refer to the former term as the "direct" term and the 
second as the "spin flip" term. We approximate V a p(j, p) 
by 



and V s { by 



with 



V a/3 (j,p) lafj U{j,p) 
V s {(j,p) « 7]U{j,p), 



(8) 
(9) 



U{j,p) 



1 



'Oj -3 



P 3 + A af3 



(1 



f - 3 



CM 



[p 2 + CM) 2 ] 3/2 



(10) 



where 7 Q ^j = d aa dpp / '(47re ), 
|/c?2;'u; 1 (2)^2(2) j d 2 2 /(47reo), di is the lattice spac- 
ing, and A a p ~ with £ the width of the Wannier 
functions, governs the short range physics. More details 
on this approximation follow shortly. In principle, we 
should let A al 3 be different for the V s { and V a p terms, 
but we will argue that we can set A a p — and maintain 
quantitative accuracy in the calculated spectra. 

We also note that the factor of F = \Jdzwl(z)w2(z)\ 
appearing in 77 is almost always nearly unity in a deep 
lattice. We denote the polarizability of state 2 by (1 +£) 
times the polarizability of state 1. In a deep lattice w a (z) 
can be approximated by a Gaussian. Then one finds F = 




For a £ = 0.4 polarization difference 



0.99. 



At long distances, these approximate interaction po- 
tentials are exact for dipoles aligned by a (perhaps very 
small) electric field, projecting the dipolar interaction 
into the lowest band, but at short distances it is just a 
convenient approximation. The form is chosen so that at 
short distances the potential goes to a constant, as does 
the true projected potential, and approximately interpo- 
lates between the exact long-distance and zero distance 
forms. In fact, we will set A a p — 0. Although for some 
calculations this can lead to divergent results, everything 
in this manuscript remains finite in the A a p — > limit, 
and comparing the A a p — > limit with alternative calcu- 
lations using finite realistic values of A a p (not presented) 
shows excellent agreement. The fundamental reason for 
this is that at low temperature and for a dilute gas, the 
physics is fairly insensitive to the short-range physics gov- 
erned by A a p. 

The p-wave interaction Hamiltonian is (the multi- 
internal state generalization of [57], as given in Ref. [4"8] ) 



1J pw — iJ pw • JJ pw ' 1J pw 



(11) 



where 



Hgft = u a pW a0 J2 [dp\(V p ^ a (i,p))4(i,p) 

-V&(*,p) (v p Vj(i,p) 

• [(y P i> p {%, p)) ip a {i, p) - p) (Vp^(», p))](12) 

with u af) = 6irti 2 b 3 a p/m, W a p = Jdz \w a (z)\ 2 \w (z)\ 2 , 
and b a /3 the p-wave scattering length between states 
a and j3. In asymptotically deep lattices the Wan- 
nier functions are approximately Gaussian, and W a p 

1- / / — / — \ ~ 1 / 2 
simplifies to W a p = ( V \T + sj W) with V °> 

the lattice depth in recoil units for state a. If the 

lattice potentials for a and (3 are the same, then 

W a8 = (l/d)v / ?i72^ 1/4 - Alternatively, the p-wave 



interaction could be written by adding a term to the 
intermolecular potential V a /3 in Eq. ^ of the form 

(2nbl / m )W p 5(p)£ s p 3 W p Jdz\w a (z)\ 2 \w f3 (z)\ 2 07j, 
where an arrow to the left indicates that the gradient 
operates to functions on the left. 
There is also a confining potential 



dpV a (i,p)ipl(i,p)i) a (i,p), (13) 



which we neglect in our first computation of the average 
spectral shifts for the homogeneous system in Sees. |III A| 
and |III B| Its effects will be included with a local density 
approximation in Sec. |III C| 

One measures the recoil-free spectrum by applying a 
probe pulse for a short time, where the probe pulse is 
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described by the Hamiltonian 



nJ2 jdpe- i ^ l (j,p}4(j,p) + H. C . (14) 



in the rotating wave approximation, with uj the probe fre- 
quency and proportional to the amplitude of the probe 
field. The spectrum is given by counting the number of 
atoms transferred to state 2 as a function of ui for a fixed 
probe time. 

As an aside, we note that in atomic clock experiments, 
the small variations of Q with the trap quantum number 
give rise to a shift even for s-wave interacting - initially 
identical - fermions [121131]. This effect is quite tiny, com- 
ing from the differences in Rabi frequencies associated 
with each trap mode: they only manifest in clocks be- 
cause of the clock experiments' extraordinary frequency 
resolution (~ 1Hz). The effect is even smaller for cur- 
rent experiments: the variation of this Rabi frequency 
for radio-frequency or microwave pulses is orders of mag- 
nitude smaller, and so the resulting shift will be much 
smaller than a Hertz. This is completely negligible in 
the current experiments. 

Finally, we comment that we calculate all frequency 
shifts relative to the level splitting of the internal states 
(e.g. rotational levels) 1 and 2 of a single molecule, S v . 
The experimentally measured spectra will then all be 
shifted by this value. 



III. SPECTRAL SHIFT: SUM-RULES 
A. General expressions for the homogeneous gas 

The spectral weight for the recoil-free spectrum is the 
number of particles transferred from the initial internal 
state in a time t, as defined above. Fermi's Golden Rule 
gives 

where the sum over i and / runs over all possible energy 
eigenstates of H, pi is the initial probability of being in 
state i, and Ej is the energy of state j for the Hamiltonian 
H. 

The average spectral frequency 



is 01 [50] 



/ dui ioIIuj) 
Jdu;I(uj) 

([H p ,H]H p ) 
<^ 2 > ' 



(16) 



(17) 



where the expectation includes both the thermal and 
quantum averages. This is an exact result (within lin- 
ear response) relating (uj) to the expectation value of an 



operator which can simply be evaluated from the defini- 
tions of Hp and H, for any Hamiltonian, and holds even 
for strongly correlated systems. 



We evaluate the commutators in Eq. (171 for each in- 
teraction term in H in turn. The spectral shift from 
s-wave interactions vanishes. 

The spectral shift from the direct and spin flip terms 
of the dipolar potential is 



j 



J dp [V 12 (j, P ) - V U (J, P) + V sf (j, p)] 
x<? 2 (j,p), (18) 



with 



(19) 



n = (V4 (0,0)^1(0,0) 
the two-dimensional density and 

g 2 (i,p) = i^(i,p)^(0 s 0)V»i(0,0)Vi(i,p)).(20) 

This holds for any two-body potential, and one sees that 
it vanishes for an s-wave contact interaction because 
(^(O)^ 1 (0)V>(0)i/>(0)) = 0. The key qualitative features 
of this result are understood by mean field theory: the 
expression, up to constant factors in each term, is the 
mean field interaction energy of the state with one molec- 
ular state changed (a superposition of states with a single 
molecular state changed at position p, summed over p, 
giving V\2 and V^f terms) minus the state with all the 
molecules in the same initial state (Vn term). 
The spectral shift for the p-wave potential is 

_ 8(u 12 Wi2 -wiiWn) 

[Upw ) — 

n 

(v^(o,p)) m(o, P ) ■ (v^(o,p))) (21) 

with ni(j, p) = tpxUi p)ipi(J> P)- The total average spec- 
tral shift is (uj) = (ojd) + (u! pw ). 



B. Weakly interacting homogeneous gas 



Equations (18) and (21) give exact, general relations 
for the spectral shift, valid in any system. Here we sim- 
plify these for weakly interacting systems. In the limit 
where the interactions induce negligible correlations, the 
system's correlations are that of a free Fermi gas. In 
other words, since we want to know the shifts to low- 
est order in the interactions, and the interaction already 
appears multiplying the operator expectation values, we 
may evaluate these expectation values to zero'th order 
in the interaction. Consequently Wick's theorem deter- 
mines the expectation values. 

The spectral width of a homogeneous gas narrows as 
the interactions decrease, as we show in Sec. |IVB| For 
experiments considered in the present manuscript, the 
spectral width is very small and the spectrum is nearly 
a perfect ^-function. In this case, the sum rule gives the 
exact location of the well-defined spectral peak. 
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1. Expressions at general temperature 

For the dipolar interaction, we first compute gi- 
Wick's theorem gives 



( G(i,p) 
\ n 



(22) 



where G(i,p) = yp\(i, p)ipi(0, 0)^}. Because there is no 
tunneling, the Green's function has the form 



with 



G(i,p) = di, G(p) (23) 
G(P) = Jdke ik "(n lk ), (24) 



where the integral runs over all the transverse momen- 
tums states of the system (which, recall, in this section 
is assumed to be homogeneous). 



The expectation values are, from Eq. (19 1, 
1 f 1 



m 
2^/3 



log (1 + e^) 



and, from Eq. (24) 
G(p) = 



_! f kdk Mkp) 



(25) 



(26) 



with Jn the zero'th Bcsscl function of the first kind. 



Equations (18), (22), and (26) give the dipolar shift 



For the p-wave interaction, Wick's theorem gives 

(Upw) = - {ui 2 W 12 - wnWn) 
n 

'(v$(0,p)) • (V^(0,p)))(n(0,p)) 

(1I12W12 - u n Wu) 
k 3 



(27T) 



X / ^ e Plk 2 /(2rn)-t,] + 1 

8m 2 T 2 



7T 



(u 1 iW 1 2-u n W 11 )U 2 (-e?i l ) 



(27) 



where Li Q is the polylog function. 



2. Thermal (Boltzmann) and degenerate gas limits 

In this section we evaluate the shifts, given above for 
a weakly interacting system at arbitrary temperature, in 
the thermal/Boltzmann (T ^ /i) and degenerate (T <C 
fi) limits. 



In the thermal limit T ^> fi the dipolar shift, given by 
Eqs. @, (pi, and ([26]) taking A aP = 0, is 



(u d ) = n (712 -Jii + Tf) 
00 

+ 2^ J dp 



dp 



1 



V//3 



1 - 3 



( mf 

\P 2 + Udi) 



3 = 1 ' 

= 2n 3 ' 2 J^n{ ll2 



(p 2 + (jdi) 2 f/ 2 
711+7?) thermal, (28) 



where n — [mj (2irj3)]e^ tl . Note the integrals in the sum 
over j vanish. The physical reason for this is that these 
integrate the electric field generated by the zero'th pan- 
cake over the plane of pancake j. Closing this planar 
surface at infinity, the total electric flux penetrating the 
surface must be zero. Since the dipolar field decays as 
1 /r 3 the contribution from the surface at infinity is neg- 
ligible, the electric flux through the plane must be zero, 
and hence the integral vanishes. 

In the degenerate limit T <C p, the dipolar shift is 



(uj d ) = n (712 - 711 + 77) 
256 



dp 



1-4 



[Jl(k FP )V 

(k FP y 



45 



nkp (712 — 711 + 77) degenerate (29) 



with kp = y/2mfj, the Fermi wavevector, and n — 
kp/(4ir) in this limit. 



The p-wave shift, Eq. (27) simplifies in the thermal 
T ^> g as to 



8mT 



u pw / 



n (U12W12 — unWu) thermal (30) 



and in the degenerate T <C p. gas to 



{UJ V 



kp 



(u 12 Wi2 — unWu) degenerate. (31) 



C. Spectra in a trap 

In this section we calculate the spectra of systems in 
a harmonic trap. We use the local density approxima- 
tion (LDA): we approximate the properties at a posi- 
tion (i, p) in the trap by those of a homogeneous sys- 
tem at an effective chemical potential p e s{j, p) = po ~ 
(muj 2 /2) [(jdi) 2 + p 2 ] with /iq the global chemical poten- 
tial of the system and to t the harmonic oscillator poten- 
tial's trapping frequency. As usual, the LDA is accurate 
when the local chemical potential variation over the cor- 
relation lengths is sufficiently small [ST]. This condition 
is met in the experiments of present interest. There is 
an additional requirement for the accuracy of the LDA, 
because we are considering a rather nonlocal quantity, 
the shifts from the the dipolar interaction energy. For 
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this quantity, accuracy of the LDA requires that the spa- 
tial variation of the local chemical potential is sufficiently 
slow compared to the length over which the dipolar in- 
teraction becomes negligible. Here, the accuracy in a 
spatially varying system will not be expected to converge 
exponentially, but rather the error e scales as a power law 
e ~ n J L (2irp)dp-^3 oc n/L in the scale L over which fi 
varies (L is on the order of the cloud width). Meanwhile, 
as seen in Eq. (18), the contributions to the shift scale 



as njl with I a characteristic microscopic scale (ther- 
mal de Broglie wavelength in the thermal gas or Fermi 
wavelength in the degenerate gas). Since (for our trap 
parameters) L is typically much larger than tens of mi- 
crometers, while I is typically a tenth of a micrometer, we 
expect this to be quite accurate. However, it will lead to 
quantitative corrections at the < 1% level. Consequently 
the LDA is a good approximation for systems with many 
particles, and is typically an excellent approximation for 
cold atomic systems. Here we take the trap to be har- 
monic and isotropic for simplicity, but these restrictions 
are straightforwardly relaxed. 

Within the LDA, the trap summed spectrum is 

T(u) = 277n 2 t^2 j ' dpn[p cS (i,p)}5 [uj - w M=Alcff(i , p) ] . 

(32) 

We can generally numerically sample the integral, ap- 
proximating the delta function as a sharply peaked 
Lorentzian of width F, 8(x) — > Sr(x) with Sr(x) = 
(l/7r)r/(x 2 + r 2 ). In some cases, however, an analytic 
formula can be derived, as we now discuss. 

The spectral weight in the frequency window (uj,uj + 
did), given by I(ui)duj, is 2nil 2 t times the number of 
atoms with spectral shifts in that window: 

l{u>) = (27rfl 2 t)4:TT[r(uj)} 2 n(r(uj))\dr/duj\ (33) 

where r(u>) is the distance to the trap center of atoms 
that have spectral shift equal to uj. Within the LDA, 
(to) as a function of r is inverted to give r(u>) and one 
similarly can calculate dui/dr. 

In the thermal gas we obtain analytic formulas. Note 
that in this regime the sole chemical potential depen- 
dence of both the p-wave and dipolar interactions comes 
from the density, which multiplies a term independent of 
the chemical potential: (uS) = Cn for some constant C . 
Moreover, in this limit, n(r) = [m/ ' {2irf3)]eP^ rAt ^ . Then 
r(w) = (l/wt)y/2/(mp)y/l3no ~ log (2npu/(Cm)) and 
\dr/du)\ = y/2mf3ojuj t f3pQ 



log (2n/3u/(Cm)) 
Thus the trap summed spectrum of the thermal gas is 



1(uj) cx W/3/io - log 



2tt/3uj 
Cm 



trap, thermal (34) 



for ui such that the term under the square root is real 
and positive, and zero otherwise. 



No simple expression exists for the case of general tem- 
perature. Nor does one exist even if the gas is deeply 
degenerate, although we mention that such formulas do 
exist for the deeply degenerate gas if either the p-wave 
or dipolar shift is included alone. For a purely dipolar 
gas in the deeply degenerate regime (T — 0), one finds 



1(oj) cx uj 1/3 



\ 



Ho I \I2-kuj 



2/3 



trap, degenerate 
(35) 



with Wdo 



^ 3 (W(7i2-7ii+r/))^ 



D. Example spectra 

Figure [3] (a) shows an example spectrum for a ther 
mal gas calculated by numerically integrating the LDA 
expression Eq. ( 32 ) with the dipolar shifts Eq 
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for a 

finite size system of 200 pancakes (although only ~ 30 of 
these are significantly occupied) and trapping frequen- 
cies and temperatures roughly corresponding to ongoing 
JILA experiments: ui t — 2ir x 100 Hz and T = 800 nK, 
with 2200 particles in the central tube. In fact, the JILA 
traps are quite anisotropic, but, since our goal here is 
not to construct detailed models of particular experi- 
ments, an isotropic trap gives an adequate caricature. 
Fig. [3] (b) shows the spectrum at T = 0, showing how 
the spectral shape changes when the gas becomes deeply 
degenerate. In contrast to the asymmetric thermal cloud 
with its small-w logarithmic divergence, a more symmet- 
ric, broad peak forms. At intermediate temperatures, the 
spectrum looks roughly like a sum of these two spectra: 
there is a broad spectral peak from the degenerate por- 
tion of the trap and a small-w logarithmic divergence. 
This shows how the spectra may be used to assess the 
degeneracy of the gas. 

The structure shown in Fig. [3] is quite general. A 
trapped system's spectra generically shows two charac- 
teristic features: a logarithmic divergence at low fre- 
quency due to thermal tails and a broad structure when 
the gas becomes strongly interacting and/or degener- 
ate [S5]. Such structure has also been seen in other sys- 
tems, such as dilute spin-polarized atomic hydrogen [T]. 

The dipolar energy scale w S caie = ^dir + w„f, with 

Wdir = (712 - TiOy 3 ^ and oj s { = V\f^0~ depends on 
the rotational states being used and the external elec- 
tric field applied. For concreteness, we estimate the shift 
in the transition from the ground rotational state |0,0) 
to the excited rotational state |1,0). The magnitudes of 
shifts for other low-lying rotational transitions are generi- 
cally the same order of magnitude. Fig.|4]plots the energy 
scale Wgcaie versus the dipole moment in units of the per- 
manent dipole moment. The permanent dipole moment 
for KRb is about 0.5 Debye. We calculate the dipolar in- 
teraction and w sca ie by computing the dipole moments of 
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FIG. 3: (a) Recoil-free spectrum with a rough modeling of 
the JILA experimental parameters. Here N e (oj)/N e (u lna , x ) is 
the number of excited particles after a probe pulse of fre- 
quency uj and duration t rescaled by the peak's maximum 
value. The overall frequency scale, related to the dipolar in- 
teraction strengths, is o; sca i c = (712 - 711 + rj) \J^0-- The 
shape of this spectrum is entirely determined by the density. 
The infinite system has a logarithmically divergent peak at 
uj — 0, but this is rounded off by a small spectral broadening 
F = 7 x 10~ 4 ixJ sca i c added to the calculation here and, less so, 
by the finiteness of the present system, assumed to be 200 lat- 
tice sites wide, (b) Recoil-free spectrum at T = (deeply de- 
generate) with LJdcg 
we have chosen fio 
plot. 



= [ m 3(10|4 (7i2 _ 7ii+r?) )2] whem 

2.6a;dcg; different fio's just rescale the 



these rotational states using a quantum rigid rotor model 
in an electric field with parameters given in Ref. [5 3) . We 
observe that this energy scale can be greater than 50 kHz. 
Then we see from Fig. [3] (a) that even for the low den- 
sities of the JILA experiments the width of the spectra 
is as large as about 5kHz. Experimentally resolving such 
scales is standard practice. 

The real part of the p-wave interaction is unknown, but 
it should be quite small at these low temperatures. On 
the other hand, the imaginary part of the p-wave scatter- 
ing, due to the reactive loss when molecules collide, can 
be much faster, and limits the cloud lifetime for KRb, 
even when it is suppressed in by using a two-dimensional 
"pancake" geometry [TS]. It gives an imaginary shift to 
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FIG. 4: The characteristic frequency of the collisional shifts 
in the thermal gas, a> sca i , for transitions between KRb's two 
lowest rotational states versus the dipole moment dn in units 
of the permanent dipole moment d p (roughly 0.5 Debye for 
KRb). Middle curve: total characteristic frequency aj sca i c = 
ui s f + oJdir- Top curve: w s f, the contribution from spin-flip 
term. Bottom curve: ojdir, the contribution from the direct 
term. 



the average frequency from the sum rule, and leads to a 
Lorentzian broadening. To determine the value of this 
broadening, note that Eq. (30 1 with 612 = is the ex- 
pectation value of the p-wave interaction energy. The 



imaginary part is the particle lifetime, since at the level 
of calculation of this equation the imaginary part of this 
energy is the quasiparticle lifetime and the only source 
of quasiparticle decay is the p-wave particle loss. Thus 
the cloud lifetime is simply this quasiparticle lifetime. 
Using cloud lifetimes measured in Ref. [TS] that decay 
is around 4/s, or a little less than a Hertz. Since the 
shift is the difference of two such rates - the 1-2 chan- 
nel and the 1-1 channel - we expect it to be at most 
on this order, and likely considerably less. Of course we 
can't rule out the possibility that the 1-2 interaction is 
much larger, but it is highly unlikely to be the orders of 
magnitude larger required to contribute substantially to 
the spectral broadening, compared to the dipolar shifts. 
Thus the p-wave broadening will be negligible under the 
current conditions of the JILA experiments. 



IV. SPECTRAL SHAPE: BROADENING 

Previously, we considered the mean spectral shifts. 
Also, because the p-wave shifts may have an imaginary 
part due to the reactive nature of certain molecules, such 
as KRb, this gives a broadening. Here we go beyond 
treating the spectra as sharp delta functions with this 
simple broadening, and consider the spectral lineshape 
for the experimentally relevant dilute gas. Broadening 
arises from several sources: differing trapping potentials 
felt by the two states due to their differing polarizabilities 
(coming from vector and tensor light shifts), collisional 
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broadening (related to the quasiparticle lifetime), tem- 
perature broadening from going beyond the pseudopoten- 
tial approximation, and higher band occupation. We will 
see that in this limit all sources of broadening arc very 
small, except the differing state polarizabilities. How- 
ever, even for current conditions, these broaden the spec- 
tra by an amount that is comparable to the width of the 
trap-summed spectra. Consequently some key features 
survive, such as the maximal shift. Moreover, by apply- 
ing the probe after aligning the dipoles at the "magic 
angle" or "magic field strength" [S3], even this broaden- 
ing can be eliminated. 

There is also the natural linewidth of the excited ro- 
tational state, given by r nat = (<5i / /c) 3 d^ 2 /(37re fi) where 
b v is the transition frequency between the excited state 
|1, 0) and the ground state |0, 0), and d\i is the dipole ma- 
trix element between these states. Even for the strongest 
dipole moments, this is ~ 1CP 10 Hz for KRb, calculating 
d\2 and 8 V in the same manner we calculated the energy 
scales above, and thus completely negligible. 



A. Differing polarizabilities 

The polarizabilities of rotational states can differ sub- 
stantially for different states due to vector and tensor 
light shifts. For example, the polarizability difference 
of KRb between the lowest and first excited rotational 
states is about 30% for lattice wavelengths giving a lattice 
spacing of 545nm [53] . In systems using optical potentials 
to confine the molecules, this leads to the molecules ex- 
periencing different external potentials. In the presently 
considered case, this means the trapping frequencies dif- 
fer by about 15% and the lattice depths differ by about 
30%. Since we are taking the lattice to be sufficiently 
deep to suppress tunneling, its only effect on physics 
along the lattice direction is to change the Wannicr func- 
tions for each state and thus renormalize the interactions. 
This was accounted for in our derivation of the lattice 
Hamiltonian. However the differing trapping frequencies 
in the transverse direction leads to broadening. 

Within the LDA, we note that if the trapping poten- 
tials differ by A(r) at location r then this shifts the spec- 
tral line of atoms at that position by s(r) — A(r) (this 
intuitive result follows formally from calculating the sum 
rule commutators of Eq. (IT)). For harmonic traps, this 



shift is s(r) = -A% + A 1 



i (a 



where u„ is the 



harmonic oscillator trapping frequency in the transverse 
direction for state a and A a are constants. Note that A a 
are positive since the optical intensity is maximal at the 
center of the trap and that there is a frequency shift there 
as well. The magnitude of the shift A\ — A^ depends on 
the exact trapping laser used, and just gives a constant 
shift to the spectrum. What is more important for our 
present concerns is the broadening coming from the sec- 
ond, spatially-dependent term. At finite temperatures, 
there is (exponentially small) occupation arbitrarily far 
from the trap center, and the shift from the second term 



increases indefinitely as one moves outwards in the trap. 
However, because the spectral density for a given value of 
this shift is proportional to density, the result of the po- 
larizability difference is a exponentially decaying tail, in 
frequency, in the spectrum. The spectrum for locations 
at distance r from the center of the trap has a polariz- 
ability difference shift 



(36) 



and the spectral weight is proportional to n(r). In the 
thermal gas, Eq. (33) with the frequency position rela- 



tion Eq. (36) thus gives the spectrum, at least in the 



tails where the polarizability difference shift is the only 
relevant one and the cloud is thermal. Thus oc 



r 2 (dr /dui)e 



yielding 



oj\e 



-0- 



er 



(37) 



where is the Heaviside step function. The spectral den- 
sity decays exponentially with frequency, with the char- 
acteristic decay frequency [(wf — )T/u:f\. Note that this 
decreases with decreasing temperature, by virtue of the 
cloud shrinking so that less particles occupy the wings, 
where the potential difference is largest. Also in the tail, 
this is the only source of shifts since densities are low, 
so this always gives the complete spectral shape of the 
spectral tail. 

For an 800nK gas with polarizabilities of the two states 
differing by 30%, this characteristic frequency scale over 
which the spectral density decays is about 4.9 kHz, com- 
parable to the spectral width of the trap-summed spec- 
tra obtained without this source of broadening. How- 
ever, we emphasize that these shifts can be eliminated 
by performing the recoil-free spectroscopy while align- 
ing the molecules at the magic angle or at the magic 
field strength, using states with nearly identical polar- 
izabilities, or using other techniques to compensate the 
polarizability difference. As such, the spectral tails can 
be a useful diagnostic of how close the trapping poten- 
tials of the two species are, and can help, for example, to 
experimentally pinpoint the magic angle. 

Figure [5] (a) illustrates the spectral shape due to polar- 
izability difference broadening, with a sufficiently dilute 
systems such that there is no interaction shift. In this 
limit, the spectrum without the polarizability difference 
broadening would simply be a (^-function. For illustra- 
tion, we assumed that uj\ < us\. For uj\ > ujf, the spec- 
trum will be reflected across the y axis. Figs. [5] (b,c) 
illustrate the effect of this broadening on the spectrum 
in the presence of dipolar interaction shifts, analogous to 
Fig. [3] 



B. Collisional lifetime 

Interparticle interactions broaden the spectral line by 
giving the quasiparticle excitations associated with the 
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FIG. 5: Polarizability difference broadening in the spectra, 
(a) Polarizability difference broadening when other shifts are 
negligible (as applies in a very dilute gas, or in the tails of the 
cloud). (b,c) Broadening with the polarizability difference 
for positive (b) and negative (c) uj\ — <jJi, with magnitude 
decreased five-fold from its true value. They compare the 
thermal spectra without broadening [blue, dashed, exactly as 
Fig. [3] (a)] with the spectra with broadening [red, solid], for 
dipole moments giving the maximal values of the collisional 
shifts relevant for the JILA parameters computed previously. 



response a finite lifetime. For weak interactions, the line- 
shape is roughly Lorentzian. We emphasize that the full 
lineshape is captured in Eq. (151. Only our subsequent 



determination of trapped system spectra by the LDA and 
treating the homogeneous systems' spectra as delta func- 
tions made an assumption that the linewidth was nar- 
row. Here we will derive the linewidth and show that 
this assumption is justified under current experimental 
conditions. Because we have found the effect of p-wave 
collisions to be substantially smaller than the dipolar col- 



lisions, we only consider the latter here. Moreover, we 
will treat only the direct term dipolar effects in a sin- 
gle pancake and restrict ourselves to the thermal gas. 
Including inter-pancake interactions and spin-flip terms 
is more tedious, and we will see that the direct, on-site 
contribution to the broadening - even at its largest ex- 
perimentally relevant values - is fairly small. Thus, we 
may expect the other effects - inter-pancake interactions 
and spin-flip terms to be small as well. We return to this 
in a bit more detail at the end of this section. 

We start by calculating the quasiparticle lifetime for 
particles in the initial rotational state. The standard 
theory of collisional lifetimes for a quasiparticle of mo- 
mentum p and frequency uj gives the decay rate |54j 

r(p,w) = J dp'dqS(uj + e p , -e p+q -e p ,_ q ) 

x [v(-q) - v(p - p' + q)} 2 
xlf(e p ')[l-f(e P+q )} [1-/(V-*)] 

- [1 - /(v)] /( e P+ g )/(V- q )l ( 38 ) 



is the Fermi function, e p 



where /(e) = ^ 
p 2 /(2m), and v(q) is the Fourier transform of the in- 
teraction potential, here the Fourier transform of Vn(p). 
This can be understood by the decay processes occuring 
in Fermi's Golden Rule, perturbing in the interactions. 
For a thermal gas, this simplifies to 



T(p,uj) = 1 J dp'dqS(ui + e p f - e 



p+q *-p"-q 



X [v(—q) - v(p - p' + q)] e 



0(V-*O 



) 

(39) 



Rather than worry about the full p and uj dependence 
of this broadening, we calculate the average broadening 
of the line and for each p treat the frequency as being 
fixed to uj = e p . With this, the typical broadening - 
that is, the broadening averaged over momenta r t y P = 
(l/iV)£ p r(p,e p )-is 



- typ 



dp 



r(p,ep). 



(40) 



Calculating r typ from Eq. (J39J) requires the Fourier 
transform of the interaction potential. The Fourier trans- 
form, neglecting An as argued to be valid earlier, is sim- 
ply 



v(q) = C - 2ir ln q + 0{q 2 



(41) 



for some constant Cq. The constant Cq cancels in 
Eq. (39). Physically, this is because it describes a 



renormalization of the contact potential, which affects 
nothing in a gas of identical fermions. Then, evalu- 
ating the integrals, the typical quasiparticle lifetime is 

r tyP = m2f ^ 3 2h " Tn. Ref. [55] argued that the recoil- 
free spectral broadening is simply the expression for the 
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quasiparticle decay rate with 711 712 — 711, giving the 
collisional spectral broadening 

T rf = —3 Tn. (42) 

It is useful to compare this broadening to the mean dipo- 
lar spectral shift: 

r m T(7i2 - 711). (43) 



47r 9/2 



Note that this ratio decreases with temperature, and is 
proportional to 712 — 711- Thus as the shifts decrease, the 
ratio of the broadening to the shift also decreases. Even 
for the largest experimentally realistic differences 712 — 
711, for example setting 711 = and taking 7 12 to be that 
associated with the permanent dipole, this ratio is about 
0.02. More typical differences will be a small fraction of 
this (< 10% or smaller). Given the already small shifts, 
this broadening is negligible, perhaps at most a hundred 
Hertz. 

So far, we have neglected the interactions with other 
pancakes and the spin flip terms. Neglecting the inter- 
pancake interaction is reasonable since the inter-pancake 
interactions are somewhat smaller than the intra-pancake 
interactions — by being further away by at least one lat- 
tice spacing. The fact that the shift vanishes for inter- 
actions between pancakes suggests that the broadening 
may be even smaller than this naive argument suggests. 
Since the effect of the intra-pancake interactions is quite 
small, it is almost certainly an excellent approximation to 
neglect the inter-pancake interaction induced collisional 
broadening. Regarding the spin flip terms, their maximal 
is less than the maximal direct term's, and thus the spin 
flip term's effect is likely to be comparably small. 



C. Broadening from interaction potential shape 

At sufficiently high temperatures, multiple scattering 
channels will become relevant, and rather than merely 
shifting the spectral line, interactions will also broaden 
it. This leads to a broadening as one leaves the cold 
collisional regime. This occurs when the thermal de 
Broglie wavelength At becomes less than or comparable 
to the range of the potential. Here, however, the tem- 
perature is sufficiently low that At is somewhat longer 
than the characteristic range of the potential, even for 
the dipolar potential. To see this, observe that presently 
At ~ 70nm while the dipole length is ~ 50nm. As ex- 
periments achieve colder temperatures, At will increase. 



D. Higher band occupation 

We have so far neglected higher band occupations. 
These have various effects, including a broadening of the 
spectra. Such broadening should mostly be, relative to 



the shifts, on the order of the relative occupation of the 
higher bands, perhaps ~ 10-20% at most here. A some- 
what larger contribution can come from s-wave collisions 
between molecules in different bands. Regardless, these 
occupations are exponentially suppressed as the tempera- 
ture is lowered, and additionally it may be possible to re- 
move them by purification techniques: e.g. higher bands 
tunnel more rapidly than the lowest band, and it may be 
possible to remove them from the lattice on a short time 
scale. 



V. 3D LATTICE 

This section considers the recoil-free spectra of 
molecules in a sufficiently deep 3D lattice such that tun- 
neling is negligible. This limit is interesting both for its 
simplicity and, as we show, its utility as a direct measure- 
ment tool for numerous important correlation functions, 
including for states created when the tunneling is signif- 
icant. Experiments on KRb molecules in 3D lattices are 
underway at JILA |56j . 

In the infinitely deep lattice, the eigenstates of the 
initial Hamiltonian, prior to application of the spectral 
probe, are \{rii}) = ®j \ n i)i where \n) i is the Fock state 
with n particles at site i. Then the total spectrum is the 
sum over sites j of delta functions whose spectral shift 
is that for exciting particles at site j; in the cigenstate 
the shift associated with exciting site j is 

K- = J^[Vi2(Rij)-V 11 (R ij ) + V a{ (R ij )]n i (44) 

where Rij is the separation of the i'th and j'th lat- 
tice sites. This is essentially the integrand/summand of 
Eq. (18), since we are considering the shift associated 



with changing the rotational state at a single site j, pro- 
jected to a 3D lattice. Each interatomic separation be- 
tween i and j will have its own unique shift, and this will 
allow us to spectroscopically measure the correlations at 
various distances. Broadening can be neglected, as the 
only source of broadening when the tunneling is zero is 
the ~Hz p-wave losses. 

Any many-body wavefunction |\&) initially all in the 
a = 1 internal state can be expanded in this basis: 



I*) 



= E A im} KM> • 



(45) 



The total spectrum is a sum of peaks with shifts given by 
Eq. (44). In particular, the total spectrum in the present 
zero-tunneling limit is I to tM = J2{m} ^{n^fnjM 
where I[ ni y(iv) is the spectrum for the eigenstate 
This spectrum is Xs ni \{oS) 
given by Eq. ( 44 ) . 



rij8(u> - with (w) 



For simplicity we consider the case where the 3D lat- 
tice is broken into 2D sheets such that the inter-layer 
interactions are negligible. This can be achieved either 
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FIG. 6: Zero tunneling recoil-free spectra, illustrating how 
to read off correlations, (a) Recoil-free spectra with near- 
est neighbor interactions only. The correspondence between 
spectral shift and spatial correlations is labeled, where the 
peak labeled (n\ n-2 713) has spectral weight that is a sum over 
all configurations and sites i with n\ particles on site i, n% in 
i's nearest neighbor sites, and 713 in the next nearest neighbor 
sites, summed over all sites i. The ". . ." indicates further cor- 
relations are irrelevant to the shift, and, independent of their 
values at longer range, all states contribute to the peak. Here, 
Zi and Z2 indicate the number of nearest and next nearest 
neighbor sites; note the maximum number of atoms on (next) 
nearest neighbor sites is Z\ (Z-z)- (b) Recoil-free spectra with 
nearest and next-nearest neighbor interactions. Longer range 
interactions will continue to split the peaks, and with increas- 
ing spectral resolution, longer range correlations can be read 
off the spectrum. Peak heights are drawn arbitrarily. 



by depopulating all but one sheet or by increasing the 
lattice spacing along one axis. The basic idea to read off 
correlations, discussed in the next paragraph, can be used 
in more general geometries, but the procedure becomes 
more complicated. 

Figure [6] shows an example of how to read correlation 
functions from the splitting of spectral peaks, explained 
in the caption. Fig. [6] (a) shows this including only near- 
est neighbor interactions, and Fig. [6] (b) includes next- 
nearest neighbor interactions. Contributions from fur- 
ther separated sites' correlations will continue the pat- 
tern of splitting the lines by smaller and smaller values. 

Thus we observe that recoil-free spectra can be used 



to cleanly extract a tremendous amount of information 
about the static correlation functions of a system in a 
3D lattice that have either zero or one particles per 
site (somewhat higher occupations can be accounted for 
straightforwardly, but lead to complications where peaks 
start to overlap). In particular, one can measure the joint 
probability distribution P(l, n\, 712, . . .), the probability 
of a site having one particle on it, ri\ atoms on all the 
nearest neighbor sites, 712 the number of atoms on all 
the next-nearest neighbor sites, and so on. Neglecting 
occupations of more than one molecule per site, the ar- 
guments above show that this function P(l, Tlx, 112, . . .) is 
simply proportional to the spectral weight of a particular 
spectral peak. 

Measuring P(l, rii, . . . , rij) requires finer spatial reso- 
lution as j increases. For a given string of occupancies 
(m, . . . ,n,j-i), Eq. (44) gives the frequency splitting be- 



tween having rij and rij + 1 particles at lattice sites that 
are j'th nearest neighbors to be 5 = V\2{i"j) — Vufa) + 
V s f(rj) where Tj is the lattice vector corresponding to 
the j'th nearest neighbor separation. The energy scales 
for KRb for current 532 nm lattice spacing square lattices 
(neglecting interactions between the 2D sheets that occur 
in a 3D lattice, and assuming that the difference in dipole 
moments for the initial state and final state can be made 
comparable to half the permanent dipole moment) enable 
one to resolve P(l, nx, n^) if one has 30Hz spectral resolu- 
tion and P(l,m,7i2,Ti3) if one can obtain 10Hz spectral 
resolution. For other molecules, the measurements can 
be even more informative. For example, with RbCs one 
can measure P(l, nx, ■ ■ ■ , n^) using a 30Hz spectral res- 
olution and P(l, ni, ■ ■ ■ > n 7) using a 10Hz spectral res- 
olution. Shorter lattice spacings enhance the range of 
correlations measurable for a given spectral resolutions. 

Note that this joint probability distribution includes 
moments, e.g. (n 2 ') and (noils), but is substantially more 
informative. In particular, one can compute any mo- 
ment between nj's that can be resolved: i.e. it gives 
(n,Q°n" :L n,2 2 ■ ■ ■ n* J } for arbitrary (atQ, ax, . . . , ctj). 

Importantly, this technique can probe correlations of 
many-body systems generated by the system with a shal- 
lower lattice, where tunneling may be important. In this 
case, one can rapidly ramp up the lattice depth from the 
value where the physics is of interest to a value where 
the tunneling is effectively zero. This projects the ini- 
tial state onto a basis of the zero-tunneling eigenstates 
|{ni}), and the dynamics during the interrogation then 
are described by the considerations above. Thus one can 
measure the joint number distribution function of the 
initial state. This is analogous to the measurement of 
on-site number statistics done in collapse-and-revival ex- 
periments [571 [SHI EH] : a shallow lattice equilibrium state 
is quenched to a deep lattice state for measurement. 

Such a technique is straightforwardly extended to mea- 
sure correlations where the initial state may have sev- 
eral rotational levels populated, and to measure inter- 
site correlation functions between atoms in these various 
states. This occurs naturally in, and would be particu- 
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larly useful for, proposals to use molecules' internal de- 
grees of freedom as "pseudospins" to emulate quantum 
magnetism [25J 03] , where the technique would probe the 
pseudospin correlations. 

Here we explicitly considered interpreting these spec- 
tra for the 2D lattice in a 2D system. This is straight- 
forwardly extended to a 3D lattice, but the form of the 
V12 — Via + V s f will be slightly more complicated due to 
interaction anisotropies. However, our discussion doesn't 
rely crucially on the form of V12 — Vi% + V s {, except for 
the fact it is long ranged. 

Finally, this measurement technique is not limited to 
molecules. For example, it may be possible to measure 
correlations of alkali atoms in optical lattices by excit- 
ing all the atoms to a Rydberg state immediately before 
the spectral interrogation. The spectral pulse would then 
chosen to drive the atoms between two Rydberg states, 
which will experience dipolar effects in a manner analo- 
gous to the molecules. 

VI. SUMMARY 

In this paper, we have calculated the recoil-free spectra 
for dipolar molecules in one dimensional lattices of two- 
dimensional pancakes, and in deep three dimensional lat- 
tices. We particularly focused on the case of driving tran- 
sitions between different rotational states. We have in- 
corporated both the dipolar interactions and the p-wave 
interactions and losses, the latter of which is important 
for reactive molecules. We gave general expressions for 
the shifts, and for the weakly interacting gas evaluated 
these expressions, as well as their simpler low tempera- 
ture (T fi, degenerate) and high temperature (T ^> /i, 
thermal) limits. We calculated the broadening of these 
lines from collisions, losses, differing polarizabilities of 
rotational states, and other mechanisms. We overviewed 
typical scales for these various effects for the case of KRb 
molecules. The overall scale for the trap summed spec- 
trum under current experimental conditions is on the or- 
der of ~ 5kHz for strong dipoles due to the dipolar in- 
teractions and the inhomogencous density of the trapped 
system. This will increase to a constant with decreas- 
ing temperature and increase indefinitely with increasing 
density (until perturbation theory breaks down). Ob- 
servation of this lineshape will confirm that the interac- 
tions in this system take the expected form and can be 
used to diagnose temperature and degeneracy. Devia- 
tions from the weakly-interacting spectral lineshapes can 
then be used to observe and characterize strong correla- 
tion physics. 

The p-wave losses induce a few Hertz broadening, and 
the real p-wave interactions should be very small. Col- 
lisional broadening is also small, probably much less 



than 100Hz for current conditions. The largest source 
of shift and broadening, other than the dipolar interac- 
tions, comes from the differing polarizabilities of rota- 
tional states. This is on the order of a few kilohertz, and 
decreases with temperature. However it can be elimi- 
nated by working, for example, at the "magic angle" or 
"magic field strength" [53]. Recoil- free spectra offer a 
useful probe of this physics. 

In the three-dimensional lattice, we showed how the 
spectra give access to a multitude of static correlation 
functions that are inaccessible by other techniques. The 
spectral weight in each peak is proportional to the prob- 
ability of a certain set of spatial configurations of a par- 
ticular set of lattice occupations. This allows one to mea- 
sure joint probability distributions for correlations among 
a few nearest neighbor sites. 

We mention that although we focused on the most im- 
mediately relevant case of dipoles aligned perpendicular 
to the pancakes and in a one dimensional lattice of two 
dimensional pancakes, most of our general formulas ex- 
tend straightforwardly to other cases, for example, three 
dimensions or other angles of alignment. Only an evalu- 
ation of the integrals changes. Presumably, many of the 
overall magnitudes typically will remain similar. 

Finally, although we focused on dipolar molecules in 
our quantitative estimates, the effects are relevant else- 
where. Rydberg atoms have a strong dipolar interac- 
tion. In alkaline earth atomic clocks, p-wave interac- 
tions and reactive losses can be important, and although 
the dipolar interactions are very small, they may be im- 
portant as clock accuracies approach their fundamental 
limits. However, we mention that we have treated the 
dipole interaction electrostatically. This is valid when the 
particle separation is much less than the wavelength of 
the photons mediating the dipole-dipole interaction. For 
molecules we expect this to be quantitatively accurate. 
However, in other contexts such as atomic clocks, the 
dynamic field becomes important when interparticle sep- 
arations are comparable to or larger than the wavelength 
of the photon associated with the dipole transition [59] . 
In this case, the treatment of dipolar interactions of the 
excited state is quite different. For example, they scale as 
1/r where r is the interparticle separation. In this con- 
text, careful choice of lattice parameters mitigates these 
effects [55]. 
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